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Abstract 

■ - - , Using the examples of pion-nucleon scattering and the nucleon mass we analyze the convergence 

^^ I of perturbative series in the framework of baryon chiral perturbation theory. For both cases we sum 

^^ ■ up sets of an infinite number of diagrams by solving equations exactly and compare the solutions 

with the perturbative contributions. 

^■ 

^ • PACS numbers: 12.38.Cy,12.39.Fe,13.75.Gx 
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I. INTRODUCTION 

Nowadays, mesonic chiral perturbation theory (ChPT) [l|, y, |3| is widely accepted as 
the low-energy theory of the strong interactions based on the underlying symmetries of 
QCD. Impressive accuracy in the description of data has been achieved in the last decade 
1^1 y; y? Ill IsL |9| ,llCi |ll|. On the other hand, the baryon sector of ChPT proved to be more 
problematic (_13]- Although there has been considerable progress in recent years (see, e.g., 
|l3J). the issue of the convergence of perturbative calculations in the nucleon sector of the 
effective theory is still of great interest. 

To compare lattice calculations with experimental data, one has to extrapolate the results 
to physical quark masses. The preferred method to cope with this problem is to calculate 
the physical quantities as functions of the quark masses in an effective field-theoretical 
approach (see, e.g., |lj, ll5|, lla, llZl)- Of course, these effective theories have a limited range 



of applicability. The aim of this work is to probe the issue of convergence of perturbative 
calculations in BChPT. To that end we consider the contribution of particular infinite sets of 
diagrams to the vrA^ scattering amplitude and to the nucleon self-energy. These contributions 
are summed up by solving integral equations analytically. For a related analysis of exact 
solutions to the Bethe-Salpeter equations in BChPT in the SU{3) sector, see Ref. [18|. 

Throughout this paper we use dimensional regularization in combination with the infrared 
renormalization scheme 19|, |20| without explicitly showing the counter-terms responsible for 



the subtractions of loop diagrams. 



II. PION-NUCLEON SCATTERING 



First, let us specify the Lagrangians required for the purposes of this work. From the 
mesonic sector we need the lowest-order Lagrangian of the SU{2) sector [2| 

£^ = —TT{d^Ud^U^) + ^^TT{U^ + U), (1) 

where f/ is a unimodular unitary matrix containing the Goldstone boson fields. In Eq. (pj, 
F denotes the pion-decay constant in the chiral limit: F^ = F[l + 0{rh)] = 92.4 MeV. Here, 
we work in the isospin-symmetric limit rriu = m^ = rh, and the lowest-order expression for 
the squared pion mass is M^ = 2Brh, where B is related to the quark condensate {qq)o in 
the chiral limit J2]. Next, the nucleon fields are collected in an isospin doublet 

n 

with two four-component Dirac fields p and n describing the proton and neutron, respec- 
tively. From the nucleon sector we need the leading-order Lagrangian omitting external 
sources 

4'i = ^ {il,D^ -m + ^iA 7m75«^) ^, (2) 

where 



and m and qa refer to the chiral limit of the physical nucleon mass and the axial-vector 
coupling constant. 

Let us consider elastic nN scattering with pi, gi the four-momenta of the incoming and 
P2, q2 the four-momenta of the outgoing nucleons and pions, respectively (see Fig. 1). The 
corresponding vertex function (amputated Green's function) can be obtained by solving the 
integral equation 

T'^P2,q2;Pi,qi) = V'''{p2,q2;Pi,qi) 

/r/"A" 
j^V'^p2,q2;p-k,k)G{p-k,k)T"'{p-k,k;p,,q,), (3) 

where p = pi + qi. V'^"' stands for the tcN effective potential and G{p — k, k) is the product 
of the (dressed) nucleon and pion propagators. Here, the effective potential is defined as the 
sum of all diagrams contributing to the vertex function, which cannot be reduced to two 
vrA^ scattering diagrams by cutting one pion line and one nucleon line. 

The standard representation for the vertex function in terms of isospin symmetric and 
antisymmetric parts reads 

p6a _ ^ha r+ + i \^^b^ ^aj p- _ ^^^ 

For our purposes it is convenient to decompose the scattering amplitude in isospin-invariant 
components 

r3/2 = r+ - r- , 

ri/2 ^ r+ + 2r-. (5) 

These vertex functions satisfy the integral equations written symbolically as 

T^ = V^ + V^GT^ , (6) 

where J = 1/2 or 3/2^ and 

y3/2 _ Y+ _Y- ^ 

^1/2 = V+ + 2V- . 

yba ^ ^ba y+ + 1 [^^^ ^a] y- _ (7) 

Suppose the potential can be written as 

v'fe.,.;p.,„)=(iA)(:;i;:|)(ij. (8) 

where the f/- depend only on p = pi + qi = p2 + q2 as is the case, e.g., for the potential^ 



-ij 



2 

V (P2,g2;Pi,gi) = -^^ (^1+^3) jp-, p2 _ ^2 • (9) 



^ Everywhere below the index / can take one of the two values 1/2 and 3/2. 

^ Eq. lO corresponds to the Weinberg- Tomozawa term plus the s-channel nucleon-pole diagram obtained 

from the Lagrangian of Eq. (0) . Note that the M-channel nucleon-pole diagram cannot be written in the 

form of Eq. ||SJ| . 
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FIG. 1: Equation for the pion-nucleon scattering amplitude. 



In this case the vertex functions T^ can also be written as 



r^(p2,g2;pi,gi) = ( i h) 



'^11 '^12 



''"21 ''"22 



Substituting Eqs. (jH)) and (fTUj) in Eq. (jH)) results in the following matrix equations 



T = V + iv gr 



where 



9 



(2 7r)" 
For the undressed propagator 

G{p - k, k) 



/il^(;)<^(-*-^-)(^n^(-S) 



m + i 0+ A;2 - M2 + z 0^ 



we obtain 



5-11 = m lNn{-P, 0) + ^ lNn{-P, 0) - /^l(-p, 0) 

^12 = ^21 = (p' + m;^) 4^1(-p, 0) - M^ /jv.(-p, 0) - Jjv : 
^22 = Mp^- ^') 4^i(-P' 0) -M^^- m) In.{-P, 0) - 



m) In ■ 



The loop integrals In-k{—P, 0), /at, and I^l{—p, 0) are given in the appendix. 
Decomposing the matrices as 

9 = 9s+i> 9v, 
and substituting the result in Eq. (|TT|l we obtain 



,2 „./ 



2 „./ 



.2 „./ 



If we define 



v,+iv^gsT^ +tp v^g^T^ +tp v^gsT^+ip v^ g^ r, , 
vl + ivl gs tI + ivlg^ tI + ivlgs t^ + ip^vl g^ r^ . 






9± = 9s± \ff'9v , 
^± — T^g 3= Y/; r^ , 



(10) 

(11) 
(12) 

(13) 



(14) 



(15) 



(16) 



(17) 
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FIG. 2: Differential cross section for 7T~p — > ir^n scattering in forward direction. The solid and 
dashed lines correspond to the non-perturbative and perturbative (tree plus one-loop order) results, 
respectively. 



and substitute in Eq. (fT^ the resulting equations reduce to the decoupled system 



Ti=vi + ivi g± tI . 



(18) 



Eqs. (Ilfij) and (fT^ are systems of matrix equations and can be solved exactly. Inserting the 
solutions in 



r^ 


= \{2r's" + rl'^) 


r: 


= ^(2rr + ry^) 


^7 


= \ [rl" - r'J') , 


-Tv 


= \ [rl'^ - rt") , 



and using the most general parity-conserving form for the on-shell T matrix, 

1 



yd 



^^ + 2 Wi+ 



B^ 



(19) 



(20) 



one can calculate the four Lorentz-invariant amplitudes as 

A^ = r^u + {M^ + 2pi- qi) t%^ + ttin r,^n + 2 (M^ + 2pi • gi) r, 

-mjv (M2 + 2pi-gi) r^^^, 
B^ = 2 r±i2 - 2 mjv t^22 + <ii - 2 m^ r^^s + {2m% + M^ + 2p^ ■ q{} r, 



'?);12 



s;;22 • 



(21) 



On the other hand by expanding Eqs. (pTjl perturbatively we can compare the results of the 
term-by-term loop expansion with the non-perturbative expression and estimate the error 
of the perturbative approximation for various kinematics. 

We calculated exactly (as closed expressions of one-loop integrals) the non-perturbative 
(re-summed) contribution and the tree plus one-loop order contributions to the n~p -^ 7r°n, 
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FIG. 3: The sum of the differential cross sections for the processes 7r~p — > Tr~p and 7r~p — > n^n in 
forward direction. The sohd and dashed Unes correspond to the non-perturbative and perturbative 
(tree plus one-loop order) results, respectively. 
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FIG. 4: Differential cross section for vr^p -^ n^p scattering in forward direction. The solid and 
dashed lines correspond to the non-perturbative and perturbative (tree plus one-loop order) results, 
respectively. 



TT p ^ TT p, and vr+p -^ iT~^p scattering processes for the potential due to the Weinberg- 
Tomozawa term, 



V'"'ip2,q2;pi,qi] 



'4F2 



L + 



(22) 



The results for differential cross sections are given in Figs. EB These figures suggest that 
the perturbative results (tree plus one-loop order) approximate the re-summed contributions 
very poorly already for s = p^ ~ m\. 
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FIG. 5: Contribution to the nucleon self-energy. 



III. NUCLEON SELF-ENERGY 



The fuU (dressed) nucleon propagator has the form 



iS{p) 



m — T, 



(23) 



where the nucleon self-energy — iE(|f) represents the one-particle-irreducible contribution to 
the two-point function. The nucleon self-energy contains the contributions of counter-terms 
so that m corresponds to the nucleon pole mass in the chiral limit. 

The physical mass m^ of the nucleon is defined as the solution to the equation 



S {rriN) = tun — m — S (m^v) = . 



(24) 



Below we calculate the contribution to the nucleon mass of an infinite number of diagrams 
shown in Fig. El The sum of these diagrams can be written in a closed form as 



-zS 



o 2 

9a 
'4^2 




dP-qi (i"g2 



(2 7r)" (2 7r)^ 



X 7^2 ■ 



- ^2- mN 
1 



r^r'"^(P2-g2,g2;Pi-gi,gi) r'^ 



- ^1- niN 



iHi 



X 



[ql -M^ + i 0+] [g2 _ M2 + i 0+] ' 
where the ttN vertex function T^"' is obtained by solving Eq. (jS)) with the potential 



V'''{p2,q2;Puqi] 



bac--c 



L + 



(25) 



(26) 



It is easily verified that E depends only on F^'^ Using the solution for F^'^ from the previous 
section in Eq. (|25|) and integrating over gi and g2 we obtain an analytic expression for the 
contribution of the diagrams in Fig. El to the nucleon mass 



5m 



o 2 

3 gA N 



4F2 D 



(27) 



where 



TV 



[TTl 



+ miy) (4 F^ — /,r) (mN — fn) I-K — {m^ — M'^ — 2m m.^ + m^) (m + m^v) Jjv^ 



D = 2mjv<8mivF^ + 4 



+ In 



[m — rriAr) /^ + [m^ — M^ — 2m m^N + m^) (m + m^v) 1^^, 



{rriN — m) It, — (m^ — M^ — 2mm]\[ + m%) (m + ttin) I^.,, 



(28) 
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FIG. 6: Contributions to the nucleon mass as functions of M. The sohd, dashed, and dashed-dotted 
hues correspond to 5mz-, Sm, and the contribution of the two-loop diagram in Fig. respectively. 



On the other hand, by expanding Eq. (|27p in powers of l/F"^ we can identify the contributions 
of each diagram separately. Using the IR renormalization scheme and substituting m = 882.8 
MeV im, m^v = 938.3 MeV, F = 92 A MeV, gA= 1.267 and M = 139.6 MeV we obtain 

Sm = -0.00233530 MeV = (-0.00230219 - 0.00003305 - 0.00000007 + ■ ■ ■) MeV . (29) 



As can be seen from Eq. ()29|) the first term in the perturbative expansion reproduces the 
non-perturbative result well and the higher-order corrections are clearly suppressed. 

It is relevant for the chiral extrapolation of lattice data to consider the nucleon mass for 
larger values of the quark masses. Therefore in Fig. IHl we plot 6m of Eq. (j28|) together with 
the leading contribution (first diagram in Fig. O and the leading non-analytic correction to 
the nucleon mass Snis = —3g\M^/{32 7iF'^) [l^] as functions of M. As can be seen from 
this figure, up to M ~ 500 MeV the non-perturbative sum of higher-order corrections is 
suppressed in comparison with the 61713 term. Also, the leading higher order contribution 
reproduces the non-perturbative result quite well. On the other hand for M > 600 MeV the 
higher-order contributions are no longer suppressed in comparison with 67713. 



IV. SUMMARY AND DISCUSSION 



In this work we have addressed the issue of convergence of perturbative calculations 
in BChPT by analyzing pion-nucleon scattering and the nucleon mass. By solving the 
equation for the nN vertex function using dimensional regularization we have obtained 
an exact expression for a sum of an infinite number of loop diagrams. The solution is 
given in a closed form in terms of one-loop integrals. We have renormalized the obtained 
non-perturbative expression by applying infrared renormalization [19|. We compared the 
perturbative contributions with the re-summed expression for elastic nN scattering. We 
find that already for s ~ m\ the perturbative results approximate the re-summed non- 
perturbative expression very poorly. As the potential, which has been iterated by solving 



8 



the equation does not receive contributions from the intermediate A state, we conclude 
that the inclusion of the A as an explicit degree of freedom does not solve the problem of 
convergence of the considered loop contributions. In our opinion, to solve this problem one 
needs to include the A degrees of freedom 22] and simultaneously consider the irN scattering 



equations. 

Next, using the non-perturbative result for the ttN scattering amplitude we have obtained 
an exact expression corresponding to a sum of an infinite number of diagrams contributing 
to the nucleon self-energy. Using the infrared renormalization scheme and comparing the 
non-perturbative contribution to the nucleon mass with contributions of the first several 
terms in its perturbative expansion we conclude that the so obtained perturbative series for 
the nucleon mass converges very well. We also considered the correction to the nucleon mass 
for larger values of quark masses and found that the re-summed higher order contributions 
become larger than the leading non-analytic contribution for M > 600 MeV. From this we 
conclude that for such values of M BChPT cannot be trusted in extrapolations of lattice 
data. Even if there are large cancelations these cannot be treated systematically in standard 
BChPT. As we have summed up only a subset of higher-order diagrams, our analysis is not 
complete and therefore our result should be considered as an estimate for an upper limit of 
the radius of convergence. 
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V. APPENDIX 

One-loop integrals: 

/j^,(-p,0) = pf^l^^{-p,0)=tj (2v,)n [(A; - py - m2 + ^0+] [F - M2 + tO+] 

[(p2 _ ^2 ^ ^2) J^^(_p^ 0)+In- I.] , 



P^ 



/Af7r(-P, 0) = i 

= 2A + 
where 



2p2 

I d^'k 1 , .^ M2 fM\ 

'- = ' ' (2vr)"P-M2 + zO+=''^''^+8^^H-J' 

d'^k 1 



(27r)'^ [{k - pf -m? + iO+][F - M^ + zO+] 



1 



167r2 



n^-m^ + M^ f M\ 2mM 
— 1^ — +^ 



-i+ " y 'M^ +^^(") 



(30) 



™n-4 C 1 1 






Vl — r22arccos(— f2), 



-1< fi < 1, 



Vfi2 - 1 In [n + v/TF^T) - ivrv/TF^I, 1 < fi, 



2mM 
Infrared renormalized expression for Int^{—p,0): 

1 p^ ~m^ + M^ 



jIR 



Nn\ 



where 

and a = — . 



-p,0) 



1 mM 



167r2 



2p2 



2 In — - 1 

\m J 



+ Fjn{n), 



(31) 



87r2 p2 

1 mA/ 

16 TT^ p2 



yr^TF 



arccos 



a+Q 



VT+2an+c^ I ' 



-1 < fi < 1, 



v/fF^I In "+^~^S^ - ^ iii|^ vrF^I, 1 < n, 



Intt — Intt{-P, 0)\p2=m% ■ 



(32) 
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